Breakdown of Hydrodynamics in the Radial Breathing Mode 
of a Strongly-Interacting Fermi Gas 
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We measure the magnetic field dependence of the frequency and damping time for the radial 
breathing mode of an optically trapped, Fermi gas of ®Li atoms near a Feshbach resonance. The mea- 
surements address the apparent discrepancy between the results of Kinast et al., [Phys. Rev. Lett. 
92, 150402 (2004)] and those of Bartenstein et al., [Phys. Rev. Lett. 92, 203201 (2004)]. Over the 
range of magnetic field from 770 G to 910 G, the measurements confirm the results of Kinast et 
al. Close to resonance, the measured frequencies are in excellent agreement with predictions for a 
unitary hydrodynamic gas. At a field of 925 G, the measured frequency begins to decrease below 
predictions. For fields near 1080 G, we observe a breakdown of hydrodynamic behavior, which 
is manifested by a sharp increase in frequency and damping rate. The observed breakdown is in 
qualitative agreement with the sharp transition observed by Bartenstein et al., at 910 G. 
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Optically trapped, strongly-interacting Fermi gases of 
atoms 0, Q are possibly the most convenient and flex- 
ible systems for exploring the rich physics of the BEC- 
BCS crossover 0, i.e., the transitional regime between 
Bose and Fermi statistics. Atomic systems are unique 
in that the state of the gas can be tuned throughout 
the whole crossover region, simply by varying a magnetic 
field. The last two years have seen explosive progress 
in experiments on degenerate, strongly-interacting Fermi 
gases which include anisotropic expansion '2j, studies of 
universal interactions 0,|l||,l5|, and molecular condensates 
on the BEG side of the crossover nil Hi]. Recently, 
microscopic evidence for superfluidity has been obtained 
by observing preformed pairs on the BCS side J, ,10] and 
by measurement of the gap [llj, 12]. The focus of this 
Rapid Communication is on extension and verification of 
experiments which probe the macroscopic properties of 
the trapped gas by studying the breathing mode [islIT^ . 
Mapping the frequency of the mode throughout the BEC- 
BCS crossover, as a function of the magnetic field, tests 
the predictions for the equation of state and can be used 
to verify current many-body calculational methods. 

Kinast et al. |l3j ]. have measured the frequencies and 
damping times of the radial breathing mode in an qpti- 
cally trapped gas of ^Li near a Feshbach resonance [Tj]. 
They find that for magnetic fields in the range of 770- 
910 G, the frequencies are close to the hydrodynamic 
predictions 15, 1 dEmilllllO] and that the damping 
rate drops rapidly as the temperature is lowered. The 
temperature dependence of the damping rate is consis- 
tent with a transition to a superfluid state at low tem- 
perature and inconsistent with a picture of a coUisional 
normal fluid [2ll 123 . |23| . Bartenstein et al., also have 
measured the magnetic field dependence of the frequen- 
cies for both the axial and radial breathing modes of an 
optically trapped gas of ^Li 0]. For the axial mode, 
the agreement with predictions for a hydrodynamic gas 
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is quite good 0, |23, 01 and the minimum damping rate 
is very low. In contrast, the measured frequencies for the 
radial mode are 9% below the predictions of hydrody- 
namics, and both the frequency of the radial breathing 
mode and the damping rate exhibit an abrupt increase at 
a magnetic field near 910 G. Further, the damping rate 
measured near 910 G by Bartenstein et al., exceeds the 
maximum damping rate allowed in a collisional gas by 
more than a factor of 5, signaling a possible transition 
between a superfluid state and a normal gas. 

The discrepancies between the two groups in the mea- 
sured frequencies and damping rates of the radial breath- 
ing mode motivated an additional study of the magnetic 
field dependence. In this paper, we describe new mea- 
surements of the magnetic field dependence of the radial 
breathing mode in an optically trapped, resonantly inter- 
acting gas of ^Li, over a magnetic field range from 750- 
1114 G (the widest range accessible to us at present). For 
fields below 950 G, we find that the measurements con- 
firm the results obtained by Kinast et al. , 0| . For fields 
near 1080 G, we observe a breakdown of hydrodynamic 
behavior, which is manifested in a sharp increase in fre- 
quency and damping rate. The observed breakdown is in 
qualitative agreement with the sharp transition observed 
by Bartenstein et al., il4], which occurred at 910 G. 

In the measurements, we prepare a highly degenerate 
50-50 mixture of the two lowest spin states of ^Li atoms 
in an ultrastable CO2 laser trap ^2^. using forced evap- 
oration near a Feshbach resonance [21 • The trap depth is 
lowered by a factor of ~ 580 over 4 s, then recompressed 
to 4.6% of the full trap depth in 1 s and held for 0.5 s to 
assure equilibrium. The number of atoms is determined 
from the column density obtained by absorption imaging 
on a two- level state-selective cycling transition 0, . In 
the measurements of the cloud column density, we take 
optical saturation into account exactly and arrange to 
have very small optical pumping out of the two-level sys- 
tem. The column density is integrated in the axial direc- 
tion and divided by the total atom number to obtain the 
one dimensional density n(a;), which is normalized to 1. 
Fitting the measured distributions of the expanded cloud 
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with a Thomas-Fermi distribution yields the cloud radius 
and the reduced temperature {T/Tp)fit ~ 0.1, where Tp 
is the Fermi temperature. 

To excite the transverse breathing mode, the trap is 
turned off abruptly (< 1 /is) and turned back on after 
a delay of = 25 fj,s. Then the sample is held for a 
variable time thoU- Finally, the trap is extinguished sud- 
denly, releasing the gas which is imaged after 1 ms. The 
temperature increase from the excitation is esitmated as 
AT/Tf < 0.015 when the gas thermalizes pj . 

To determine the breathing mode frequency and damp- 
ing time, the transverse spatial distribution of the cloud 
is fit with a one dimensional zero-temperature Thomas- 
Fermi (T-F) profile to extract the radius, axF as a 
function of thoid- To acquire each curve, 60-90 equally 
spaced values of thoid are chosen in the time range of 
interest. The chosen values of thoU are randomly or- 
dered during data acquisition to avoid systematic er- 
ror. Three full sequences are obtained and averaged. 
The averaged data is fit with a damped sinusoid xq + 
Aeyi-p{-thoid/T) sin{u! thoid + For a to = 25 ^s exci- 
tation time, we observe damping times up to r = 7 ms 
for both the interacting and noninteracting gas at 4.6% 
of full trap depth. 

We normalize the breathing mode frequencies to the 
transverse oscillation frequency uj± for the noninteract- 
ing gas at 526 G, where the scattering length is nearly 
zero |23| ■ Parametric resonance measurements yield 
the trap oscillation frequencies (uncorrected for anhar- 
monicity) oj^ ^ 2-^ x 1600(10) Hz, ujy ^ 2-^ x 1500(10) 
Hz and = 27r x 70(5) Hz at 4.6% of full trap depth. 
As in our previous study, the frequency uj^ is also cal- 
ibrated by exciting the breathing mode in a noninter- 
acting sample, which yields agreement to < 1%. In ad- 
dition, the measured breathing mode frequency 2 uox is 
used to scale the value of lo±_ for measurements made on 
different days, since small changes in the laser power al- 
ter the trap oscillation frequencies. Typically, we obtain 
a total atom number A'^ = 2.0(0.2) x 10^ at tempera- 
tures {T/Tp)fit ~ 0.1. From the measured trap frequen- 
cies, we find (after correction for anharmonicity using 
Eq. n below) u)±_ = ^ujxOJy = 27r x 1596(7) Hz, and 

uj — {ojxUjyUJzY^^ = 2tt X 560(13) Hz. For these parame- 
ters, the typical Fermi temperature for a noninteracting 
gas is Tp = {SNy^'^hid/kB — 2.2 ^K, small compared to 
the final trap depth of 35 /iK. 

Anharmonicity arising from the gaussian trapping po- 
tential reduces the measured frequencies from their 
zero-energy harmonic oscillator values uj. The transverse 
oscillation frequency uj±_ for the noninteracting gas at 526 
G is corrected using the formula p^ . 
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w here U is th e trap depth, M is the ^Li mass, and bs — 
■\/l + {ujxty = 10-3 is the ballistic expansion factor for 
the release time t = 1 ms. Note that the rms width 
for the ballistically expanding gas is XBrms = cr^p/y/S, 



where a^p is the transverse T-F radius of the ballistically 
expanded cloud at 526 G. 

The measured hydrodynamic frequencies ujh are cor- 
rected for anharmonicity using |28| 



loh/uT^^ = 1 + (32/25)yi073Mc.ix?„J(6l,C/). (2) 

Here, bH is the hydrodynamic expansion factor |2^. 
11.3 after 1 ms, and Xrms = cttf/VSi is the rms width 
for the interacting gas after expansion. Note that bn 
rather than 6_b is used to estimate the trapped cloud rms 
radius because we observe anisotropic expansion of the 
gas ,2] which is in good agreement with the predictions 
of hydrodynamics [3, ll^ ■ 

The extracted ratios of the corrected frequencies, 
^h/<^±, and the measured damping ratios 1/{uj±t) are 
given in Table 1 as a function of magnetic field. The 
frequencies are plotted in Fig. as a function of mag- 
netic field. The lower horizontal scale gives l/{kpa), 
where kp is the Fermi wavevector at the trap center 
hkp = \J2M ksTp and a = a{B) is the scattering length. 
Note that we assume that the Feshbach resonance is 
located at 834 G, consistent with the best recent esti- 
mates ||32i,i3l|. The solid curves in Fig.^show the predic- 
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FIG. 1: Magnetic field dependence of the frequency tJ of the 
radial breathing mode. Solid line is the theory based on su- 
perfluid hydrodynamics from Hu et al. p^ . Dashed line is 
the hydrodynamic frequency for a unitarity gas, \J 10/3, as 
predicted at resonance. Note that the top (magnetic field) 
axis is not linear. 

tions for the frequency by Hu et al. based on assumptions 
of superfluid hydrodynamics 0, • 

The data exhibit several interesting features. For fields 
near resonance, the measured radial breathing mode fre- 
quencies are in verygood agreement with hydrodynamic 
theory pjl ITtI fla . [l9j | , confirming our previous mea- 
surements fl3|. Further, the measured breathing mode 
frequency at 840 G, loh/^jJl = 1.829(.006), is very close 
to that predicted for a unitary-limited, hydrodynamic 
Fermi gas, where ljh/^i^ — \/l0/3 = 1.826. Fig. |3 
shows that the damping rates are small near resonance. 
The maximum damping time of 7 ms, obtained near 
resonance, corresponds to a minimum damping ratio 
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1/{lo±t) = 0.014, or about 20 periods of hydrodynamic 
oscillation. 
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FIG. 2: Magnetic field dependence of the damping rate 1/r of 
the radial breathing mode. Damping rate (solid circles, axis 
labelled on the right); frequency (data-open circles, theory- 
solid curve, axis labelled on the left). Solid curve (predicted 
frequency). The maximum damping ratio occurs at 1080G 
and exceeds the maximum 0.09 permitted for collisional dy- 
namics by a factor of two. Note that the top (magnetic field) 
axis is not linear. 

Away from resonance, the data in Fig. ^ exhibit qual- 
itative agreement with the hydrodynamic theory of Hu 
et al. [131 . However, the frequencies measured well be- 
low resonance, at 750 G and 770 G, are lower than pre- 
dicted by approximately two standard deviations. Fur- 
ther, above resonance, the measured magnetic field de- 
pendence of the frequency is compressed. In addition, 
there is a region where the data are consistently below 
the theory and this is followed by a region where the data 
are consistently above. The frequency reaches its lowest 
value at 925 G, dropping below the theory by three stan- 
dard deviations, and then it rises rapidly at 1080 G. At 
the highest field achievable with our magnets, 1114 G, 
the frequency is well above the theoretical prediction, 
but does not reach the asymptote, 2lOx, for a noninter- 
acting gas. In this region above resonance, we observe 
a clear breakdown of hydrodynamics, both in the rapid 
change in frequency at 1080 G and in the damping rate. 
A transition is obvious in Fig. |2 which shows the damp- 
ing rate as a function of B. Near 1080 G, the damp- 
ing rate increases to 0.17 wj^, well above the maximum 
0.09 attainable by collisional hydrodynamics [2lll2^. 
The observed behavior is therefore in qualitative agree- 
ment with the abrupt frequency change and increased 
damping observed by the Innsbruck group [T^ . However, 
our results differ from those of Bartenstein et al. , in two 
important respects: i) Near resonance, our observations 
for the radial breathing mode agree with hydrodynamic 
theory; ii) the field of 1080 G, at which we observe the 
complete breakdown of hydrodynamics and rapid damp- 
ing, is substantially higher than the 910 G value observed 
by the Innsbruck group p^ . 

The Innsbruck group has suggested that the break- 



down may arise when the zero temperature BCS energy 
gap A is smaller than the collective mode energy 
Zubarev 0, has determined the trap averaged value 
of A and shows that if the abrupt transition occurs at 
910 G for the Innsbruck group, then for the conditions 
of our trap, we should observe the transition near 1000 
G, assuming 3 x 10^ atoms. Using similar estimates, we 
find that at 910 G, A ~ 2?iw for the conditions of the 
Innsbruck group and A ~ lioj for our trap conditions at 
1080 G. Since single particle excitations have a minimum 
energy of 2 A, perhaps the condition A = Tiuj underesti- 
mates the required collective mode energy for complete 
pair breaking. However, it has been noted by Heisel- 
berg [s^ that in-gap surface modes have a smaller energy 
than A, but scale similarly, possibly accounting for the 
observed breakdown at A > Tilo. 

For the conditions of the Innsbruck group, the break- 
point occurs at the same magnetic field of 910 G at a shal- 
low t rap depth and for a trap depth increased by a factor 
of 9 hJ. This magnetic field independence was first ex- 
plained by the Innsbruck group and arises as follows. In 
their case (but not in ours) the axial confinement is pro- 
vided primarily by the bias field magnets, and the ratio 
of the Fermi energy to ?iu;j_ and hence to fioj decreases 
with increasing trap depth. Then, using the most recent 
values of the scattering length 6Y\ , we find the increase 
in fci? in the exponent of the pairing gap compensates for 
this reduction at 930 G, consistent with their argument. 

We can obtain an alternative estimate of the magnetic 
field at which an a brup t change occurs using an idea due 
to Falco and Stoof jSJ • They suggest that when the Fes- 
hbach {v = 38) molecular state in the singlet potential 
has a higher energy than the two-particle Fermi energy 
(relative to the zero of the triplet potential), the system 
becomes BCS-like. In this case, the system may not be 
interacting strongly enough to remain superfluid at the 
temperatures we achieve. A simple estimate of the break- 
down magnetic field is then obtained from 



2(6p(x)). 



(3) 



We assume that near resonance, interactions make the 
same contribution to the chemical potential of free triplet 
atoms and to very weakly bound atoms in large singlet 
molecules. We also assume the energy of the molecular 
state is eb = fi^ /Ma? relative to the zero of energy in 
the triplet potential, which is valid in the two-body case 
near the Feshbach resonance. Using the trap averaged 
local Fermi energy, we obtain 2(eF(x)) = (5/4) ksTp = 
{5/8)??Pp/M, so that l/ikpa) = -^/b/S = -0.79. 
From Tabled we see that at B ^ 1080 G, l/ikpa) = 
—0.74, in reasonable agreement. Unfortunately, our sim- 
ple hypothesis cannot explain the insensitivity of the B 
field value with respect to trap depth observed by Barten- 
stein et al.: Scaling our Fermi energy down from 2.2 /xK 
to 1.2 /iK for the conditions of the Innsbruck group, we 
predict a transition near B = 970 G, higher than ob- 
served. 
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B(G) 


l/ikpa) 




l/(c^xr) 


750 


0.728 


1.904(.015) 


0.095(.009) 


770 


0.504 


1.887(.005) 


0.071(.007) 


815 


0.124 


1.855(.017) 




840 


-0.034 


1.829(.006) 


0.014(.001) 


850 


-0.090 


1.816(.014) 


0.016(.005) 


860 


-0.142 


1.817(.015) 




870 


-0.190 


1.815(.012) 


0.017(.004) 


880 


-0.235 


1.807(.015) 




910 


-0.354 


1.779(.015) 




925 


-0.405 


1.759(.008) 


0.034(.006) 


969 


-0.531 


1.766(.007) 


0.058(.007) 


975 


-0.546 


1.778(.007) 


0.059(.007) 


1000 


-0.603 


1.763(.018) 


0.087(.016) 


1025 


-0.652 


1.786(.011) 


0.071(.010) 


1035 


-0.670 


1.781(.015) 


0.091(.016) 


1050 


-0.695 


1.764(.033) 


0.098(.020) 


1065 


-0.719 


1.786(.016) 


0.142(.020) 


1080 


-0.740 


1.796(.021) 


0.166(.028) 


1090 


-0.754 


1.803(.021) 


0.126(.016) 


1114 


-0.783 


1.823(.024) 


0.083(.021) 



TABLE I: Breathing mode frequencies u and damping rates 
1/r as a function of applied magnetic field B. Error estimates 
include the statistical error from the fit only. * Denotes data 
taken with a 50 /iS excitation time [l^ for which we omit the 
damping rate. a{B) is determined from Ref. [3ll | 



In summary, we have measured the magnetic field de- 
pendence of the frequencies and damping rates of the 
radial breathing mode for a strongly interacting Fermi 
gas of ^Li. The measurements are in very good agree- 
ment with the predictions of hydrodynamics near the 
Feshbach resonance at 834 G, but we observe a break- 
down of hydrodynamics away from resonance. Beginning 
at higher fields near 925 G we observe first a frequency 
decrease, well below predictions, and then an abrupt fre- 
quency increase at a field of 1080 G, accompanied by 
rapid damping, which exceeds the maximum damping 
for a coUisional normal fluid. 
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